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Abstract 



We propose a model for the description of the transport properties of metallic films 
on a large scale of slab thickness. This model is based on solving the linearized 
Boltzmann equation in relaxation-time approximation using ab initio calculations 
within the framework of the density functional theory. The expression for the re- 
laxation time is derived from the microscopic treatment of the scattering processes 
and provides the correct thickness dependence for very thin as well as very thick 
films. The method is applied to the calculation of the in-plane conductivity and the 
Drude-type plasma frequency of thin Cu(OOl) films in the thickness range between 
1 and 32 monolayers. 
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The electronic structure of ultrathin metallic films is of great scientific interest. 
Their transport properties can substantially deviate from the bulk behaviour 
owing to classical and quantum size effects and the existence of new types of 
electronic states (e.g. surface states). Therefore, the study of these systems is 
very important for the microscopic understanding of transport in nanoscaled 
electronic devices. 



* Corresponding author 

Email address: zahn@physik.uni-halle.de (Peter Zahn). 

Preprint submitted to Elsevier Science 2 February 2008 



Classical size effects in the conductivity as explained by the Fuchs-Sondheimer 
theory occur if the slab thickness d is comparable to the bulk mean free path 
[1,2]. In the limit of small d compared to the mean free path this model predicts 
a linear thickness dependence of the slab conductivity. The application of the 
model to this case is quite questionable, since a three-dimensional electron 
motion is assumed. Nevertheless, this behaviour is confirmed by ultrathin film 
experiments. 

When the thickness d approaches the Fermi wavelength A^ quantum size ef- 
fects (QSE) occur which cause oscillations in typical features of the electronic 
structure, e.g. density of states at the Fermi level, as a function of slab thick- 
ness. Regarding the transport properties Trivedi and Ashcroft predicted a 
saw-tooth like oscillation of the in-plane conductivity with thickness [3]. In 
their free-electron model the oscillation period is equal to half Xp. In real- 
ity, the situation is more complicated because of the discrete increase of the 
film thickness as well as the change in the electronic structure during growth. 
There are several interesting experiments which demonstrate QSE. For exam- 
ple, the observation of QSE in the dependence of the Hall coefficient and the 
specific conductivity on the thickness of ultrathin Pb(lll) films [4], the step 
height oscillations during layer-by- layer growth of Pb on Ge(OOl) [5], and the 
manifestation of quantum-well states in angle-resolved photoemission from Ag 
films [6] are worthy of mention. 

Recently, the conductivity of ultrathin Pb(lll) films consisting of N = 1 . . .8 
monolayers (ML) was studied by means of ab initio calculations using the full- 
potential linearized augmented plane- wave method [7]. The results obtained 
are in a good agreement with the experiment [8], and showed a correlation 
between QSE oscillation in the conductivity and in the density of states at 
the Fermi level. 

In an earlier paper, we have investigated the evolution of the Co/Cu multilayer 
conductivity during growth using the Boltzmann equation in relaxation-time 
approximation [9]. The transport properties of Cu(OOl) films of 2 to 5 free- 
standing ML have already been studied. Now this technique is extended to Cu 
slabs in the thickness range between 1 and 32 ML. We propose a model for 
the relaxation time of electrons caused by elastic impurity scattering in slabs 
which allows us to describe the conductivity in the limit of both very thin and 
very thick films. 



2 The model 



The electronic structure of the considered systems is calculated self-consistently 
using a screened Korringa-Kohn-Rostoker (SKKR) method [10,11,12,13,14]. 



For our calculations a screening potential with a barrier height of 4 Ry is 
used and the screened structure constants include coupling to the next four 
nearest neighbours. Spherical potentials in the atomic sphere approximation 
(ASA) are used. The calculations employ the exchange-correlation potential 
in the LDA of Vosko, Wilk, and Nusair [15]. An angular momentum cut-off 
at Imax = 3 is used for the Green's function, thus implying a cut-off for the 
charge- density components at 2l max = 6. Due to the slab geometry the states 
at the Fermi level are described by lines in the two-dimensional (2D) Brillouin 
zone (BZ). The Fermi surface integrals have to be replaced by line integrals 
[9]. In the considered geometry the irreducible part equals one eighth of the 
2D BZ and contains in our calculations about 28000 k||-points to reach conver- 
gence. For comparison, the Fermi surface integration for the bulk Cu system 
is performed on about 350000 k-points in the irreducible part of the fee BZ. 
We neglect the lattice relaxations at the surfaces. The fee lattice constant is 
taken equal to A=6.76 a.u. [16]. 

The transport is treated quasiclassically by solving the linearized Boltzmann 
equation which in the anisotropic relaxation-time approximation gives the 
following expression for the components of the conductivity tensor a of a 
crystal [17,18] 
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where e is the electron charge, V is the system volume, Ep the Fermi energy, 
v k the component of the group velocity with a Cartesian coordinate % of an 
electron in the band n with crystal momentum vector k and relaxation time 
Tfc at the energy E k (k is a shorthand notation for k and n). The factor 2 is 
due to the spin degeneracy. For a cubic crystal and a 2D lattice with square 
symmetry the transport properties are isotropic a^ = crSij. 

In the case of film geometry the relaxation time Tfu m has to be determined by 
the scattering at impurities inside of the slab as well as on surface imperfections 
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Since our calculations are performed at T = OK, the inelastic scattering by 
phonons is neglected. First, we consider the scattering rate (t 6 ) -1 caused by 
bulk-like defects inside the slab. The microscopic transition probability for an 
electron to be scattered from a state k into a state k', is given by Fermi's 
Golden Rule [10,18] 

2tt 
Pkk' = -^r\Tkk'\ ${E k — E k i) , (3) 
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where the transition matrix T kk < is defined by 

T kk ,=<4/ k \AV\^ k ,> . (4) 



Here AV denotes the perturbation of the potential, ^^ and ^k' are the un- 
perturbed and perturbed Bloch states of the system, respectively. Within the 
assumption of 5-scatterers [10] with the scattering strength tb at all impurity 
sites rj in the slab 

M b 

AV(r) = 5>*(r - r 4 ) , (5) 

the relaxation time r\ is given by the expression [10] 

1 2vr Mb o 

- = J2 P kk' = ^t 2 b J2\^k(r l )\ 2 n(r l ,E F ) . (6) 

T k k' n i=l 

The number Mj, denotes the number of bulk-like defects in the slab. Equation 
(6) includes the probability amplitude of a state k at the impurity position 
|\&fc(rj)| 2 and the local density of states at Fermi level n(r^, Ep). Further sim- 
plifications can be achieved by averaging the scattering rate (6) over all states 
k at the Fermi energy Ep and over the configurations of the impurity sites 
IV We assume the same probability amplitude of perturbed Bloch state at 
all impurity sites in the film. Introducing the bulk density Qb = Mb/ Sd of 
scattering centres, with S the area of the film and d the slab thickness, the 
averaged relaxation time is given by 

1_ _ 2-Kn f u m {E F )tlQb , . 
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where nfu m (Ep) is the total film density of states at the Fermi level (per slab 
unit cell) . The area of the film unit cell is given by the film area divided by 
the number of unit cells in the slab So = S/N ce u. 

Now, to describe the second term in Eq.(2) we introduce the additional number 
M s of the surface imperfections on the slab surfaces with the surface density 
g s = M s /S. Similar to the above derivation, an expression for t s is obtained 

1 27in surf (Ep)t 2 s g s 
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Here n sur f(Ep) is the local density of states in the surface region at the Fermi 
level (per slab unit cell), t s is the parameter for the scattering strength of the 
surface imperfections. D s is the thickness of the surface region which contains 
the surface defects. 

Introducing the defects of scattering strength tj, in a bulk material, the bulk 
scattering rate is obtained similar 

1 _ 27mb u ik(E F )t 2 b Qb , v 
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where ribuik(E F ) is the density of states at the Fermi level per bulk unit cell. 
The volume of the bulk unit cell is expressed by SqD with the interlayer 
distance D = d/N, where N denotes the number of monolayers in the slab. 

Now, introducing the dimensionless function 

7(A n = T jJL= t 2 sQsn surf (E F )N ^ 
t s t 2 b Q b n film (E F )D s 

and using equations (7), (8), and (9), the relaxation time in the slab is obtained 
by 

__ Nnbuik(E F ) T bu i k 

TfUm ~ n film (E F ) ' (1 + 7 (A0/A0 • [ii) 

The function 7 (AT) has a clear physical meaning. It defines the ratio between 
the strength of the electron scattering at surface imperfections and at impuri- 
ties inside of the slab. By virtue of the proportionality of nfu m (E F ) to A" the 
function 'y(N) is approximated in this work by a constant value 7, even for the 
case of small thicknesses. Moderate changes in 'y(N) for small thicknesses may 
change the results quantitatively, but the general trend remains unchanged. In 
the limit of large A" ^> 7 the relaxation time in the film Tfu m reaches the bulk 
value Tbuik- In the limit of very thin films the expression (11) is equivalent to 
the one used in the ref. [7] and that derived from a free-electron model assum- 
ing elastic and diffuse surface scattering [19,20]. The authors of [7] considered 
ultrathin Pb(lll) films (N = 1 . . .8 ML) with only one adjustable parameter 
like t 2 s g s in Eq.(8), and they obtained a linear dependence of the conductivity 
on the slab thickness which is in a good agreement with the experiment [8] for 
the considered thickness range. In addition, we have introduced t 2 Qb to treat 
thicker slabs. Their conductivity is appreciably determined by the properties 
of the bulk region of the slab and does not exhibit the linear thickness de- 
pendence (see, e.g. [4]). Normalizing the film conductivity (Tfum to the bulk 
conductivity <?buik it is possible to keep one adjustable parameter 7 (AT) only. 

Finally, performing the integration in Eq.(l) over the Fermi surface (for the 
slab — the Fermi lines) [21], we obtain in case of cubic symmetry 

(^ film _ 37rn bu i k (E F ) J2 n I dl n \v n (E F )\ 
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where dl n and dS n are differential Fermi line and Fermi surface elements, 
respectively. 

In case of constant relaxation times Tfu m = Tbuik f° r & U slab thicknesses (like 
in ref. [9]) the ratio of the conductivities <Jfu m /<Jbuik can be transformed to the 



ratio w-pjiim/wp^buik 

It, ,,,... / StV f^I l„ CRn^l \ ^ 

(13) 



Up, film __ ( 37rJ2nJdln\Vn(E F )\ V 2 

up, bulk \NDJ2nJdS n \v n {E F )\J 



where wpjum an d ojp,buik are the slab and the bulk Drude-type plasma frequen- 
cies [22], respectively. This ratio can be extracted from infrared transmission 
experiments [23]. 



3 Results 



In Fig.l we show the total slab density of states at the Fermi level nfu m (EF) 
determined from the one-particle energy spectrum E n (k). In addition to the 
peak for films of 2 and 3 ML in thickness large densities at the Fermi level are 
obtained for the thicknesses of 9, 17, and 30 ML. In the asymptotic case of 
very large N the oscillations in ti(Ef) are governed by Fermi surface nesting 
vectors along the (001) direction, which are known to dominate the oscilla- 
tions of interlayer exchange coupling [24]. In the investigated thickness range 
the oscillations with a period of 2.6 ML could be identified by Fourier anal- 
ysis. This oscillation corresponds to the short-period oscillation of interlayer 
exchange coupling [24]. 

The evolution of electrical conductivity with film thickness, Eq.(12), is pre- 
sented in Fig. 2 with different values of the parameter 7 for the strength of 
the surface scattering. In principle, it is possible to get the value of 7 (or the 
functional dependence 'j(N)) from a comparison with experimental data. Un- 
fortunately, experimental data concerning the conductivity during growth of 
very thin Cu slabs are missing. The reason is the prefered initial growth mode 
of Cu on insulating substrates, which forms islands instead of a continous 
film [25]. Reasonable values of 7 can roughly be estimated from a comparison 
with the experimental data for the conductivity of Fe [22] and Pb [4] films. 
The ratio <Jfu m /<Jbuik is about 0.6 for an Fe(001) slab with 31 ML (f» 4.5 nm) 
thickness, and <7fu m /<Tbuik is about 0.3 for a Pb(lll) film with 18 ML (« 5.1 
nm) thickness, assuming a bulk Pb conductivity of 0.213 (/ificm)" 1 [26]. Of 
course, these data show the strong dependence on the material. Nevertheless, 
the range of 7 should be between 10 and 100 to describe the existing experi- 
mental situations. In this range, strong QSE's can be recognized from Fig.2, 
which are clearly correlated to the behaviour of n(Ep) in Fig.l. 

Figure 3 shows the results for the ratio of plasma frequencies 0Jpju m /ujp,buik- 
From this picture a 1/JV behaviour for large N is evident. In the framework of a 
model of free electrons in a box with the height d = ND and Fermi wave vector 
kp the asymptotic behaviour for large N is obtained as {ojpjum/up 



\2 r^ 

bulk) ~ 



1 — 37r/4kpd — 7r 2 /4k F d 2 . The same expression was found for the asymptotic 
behaviour of the charge density in ref.[3], since up is proportional to the charge 
density in the free-electron model. 

The strong increase of wpju m for small N is related to the increase of the 
Fermi velocity caused by the large contribution of fast surface states. This is 
demonstarted in Fig. 4 presenting the Fermi velocity vf = (\vk{EF)\)k aver- 
aged over all states k. The ratio of vp for a 32 ML Cu(OOl) slab to the averaged 
Fermi velocity of Cu bulk is about 0.766. It is obvious, that this ratio does 
not reach one even for large N, since the Fermi velocity in the slab geometry 
corresponds to the projection of the bulk Fermi velocity into a plane perpen- 
dicular to the slab growth direction. The free-electron model gives a ratio of 
tt/4 ps 0.785. The strong similarities of our results calculated for Cu slabs 
with those found from a free-electron model are direct consequences of the 
free-electron character of the Cu Fermi surface [18]. 

The qualitative behaviour of the ujpju m /ujp t buik is very close to data obtained 
from infrared transmission spectroscopy performed during the evaporation of 
Fe on MgO(OOl) in ultrahigh vacuum (see fig. 5 in ref. [23]). Although we have 
to note, that in the experiment the strong increase of cjpju m for very thin slabs 
can be due to a strong influence from the morphology and the large roughness 
of the films just above the percolation threshold. We should mention also, that 
up, film does not saturate to a constant value after 4-5 ML as we reported in 
an earlier paper [9]. This conclusion was drawn based on the results of Cu 
slabs with thickness of from 2 to 5 free-standing monolayers. 



4 Conclusion 



We present a model for the description of the transport properties of thin 
metallic films based on the solution of the linearized Boltzmann equation with 
a quite simple and at the same time reasonable assumption for the relaxation 
time. This approximation for the relaxation time derived from the microscopic 
treatment of the scattering processes has only one adjustable parameter, and 
is valid in the limits of very thin as well as very thick films. The qualitative 
behaviour of the slab conductivity found in this approach is close to the ex- 
perimental data for thin metallic films. For very thin films we see a linear 
dependence of the conductivity on the thickness which was expected from the 
Fuchs-Sondheimer theory. In addition, the influence of QSE's is superimposed 
and causes the dips on the otherwise monotonic curve. The calculated Drude- 
type plasma frequency shows the dominance of the surface electronic structure 
for thicknesses smaller than 5 ML . So, both classical and quantum size effects, 
as well as surface states manifest in our results. 



The method applied in this paper can be used for metallic films with arbitrary 
thickness limited by the computer power only. 
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Fig. 1. Density of states at the Fermi level per atom as a function of Cu layer 
thickness. 
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Fig. 2. Film conductivity (in units of the bulk conductivity) as a function of Cu 
layer thickness. 
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Fig. 3. Film plasma frequency (in units of the bulk plasma frequency) as a function 
of Cu layer thickness. 
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Fig. 4. Averaged Fermi velocity of the Cu slabs as a function of thickness. The 
dashed line marks the bulk value. 
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